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ABSTRACT 

We study the response of giant stars to mass loss. One-dimensional simulations of red and asymp- 
totic giant branch stars with mass loss rates from 10 -3 up to a few M Q /yr show in no case any 
significant radius increase. The largest radius increase of 0.2% was found in the case with the lowest 
mass loss rate. For dynamical-timescale mass loss rates, that may be encountered during a common 
envelope phase, the evolution is not adiabatic. The superadiabatic outer layer of the giant's enve- 
lope has a local thermal timescale comparable to the dynamical timescale. Therefore, this layer has 
enough time to readjust thermally. Moreover, the giant star is driven out of hydrostatic equilibrium 
and evolves dynamically. In these cases no increase of the stellar radius with respect to its initial 
value is found. If the mass loss rate is high enough, the superadiabaticity of the outer layer is lost 
progressively and a radiative zone forms due to a combination of thermal and dynamical readjustment. 
Conditions for unstable mass transfer based on adiabatic mass loss models that predict a significant 
radius increase, may need to be re-evaluated. 

Subject headings: binaries: close — binaries: general — methods: numerical — stars: evolution — 
stars: general — stars: mass loss 



1. INTRODUCTION 

Understanding how stars respond when they lose mass 
is a key ingredient on which binary evolution models de- 
pend. This response is particularly important in the con- 
text of interacting binaries with a donor filling its Roche 
lobe on a giant branch. In the oversimplifying case of 
conservative mass transfer, the orbital separation shrinks 
if the giant donor is more massive than the companion. If 
in the meantime the giant star expands or does not con- 
tract faster than the orbit shrinks, this positive feedback 
leads to an increase of the mass transfer rate. The stellar 
response to mass loss therefore dictates, along with how 
angular momentum is lost, wheth er or not a given system 



enters a common envelope phase (Paczynski|1976 ). Con- 



a gr 



sequently, it significantly af fects the re sults of population 
synthesis studies (see, e.g., Politano et al.|[2010 ). 

Also, detailed 3D hydrodynamical simulations of the 
dynamical common envelope phase have shown that a 
giant's simulated envelope material is significantly lifted, 
but most of it does not reach escape velocity under the 
present model i ng assumptions (Pas sy et al.|2012 Ricker 
fc Taam|2012| . |De Marco et al. ( |201l[ ) suggest ed that an 



expansion of the g i ant as a result o f mass loss (Hjellming 
fc Webbink] |1987j |Ge~et al.||2010[ ) might contribute to 



the envelope ejection. Such an expansion of mass-losing 
giants was rec ently questioned in a Letter by [Woods Ez\ 
Ivanova| ( |2l)iT] ). 

T herefore, we study the radius response of mass-losing 
giants again, with detailed microphysics, using the one- 
dimensional stellar evolution code MESA (Module for 



Experiment in Stellar Astrophysics, |Paxton et aT]|2011 ) 
Such a tool - although it neglects three-dimensional ef- 
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fects - allows to remove one or more simplifying assump- 
tions adopted in some previous studies, which we men- 
tion in the next paragraphs. 

For a star following a polytropic stratification of in- 
dex n with an adiabatic index 7 = 1 + 1/n, the Lane- 
Emden equation leads to the mass-radius relation be- 
tween a standard solution of radius Rq and mass Mq, 
and a perturbed polytrope of radius R and mass M: 



R 
Rq 



M 
Mo 



(1) 



Note that Equation Q is only valid for a fixed adiabat 
throughout the stellar interior. For the comp lete deriva- 
tion, see, e.g., |Hjellming fc Webbink (1987) or Carroll 
& Ostlie (|2006p. For an ideal gas equation of state, the 



specific entropy follows a simple expression: 

s(m) = s + (1 + 1/n - i)c v ln(/j) (2) 

where sq is a constant and c v is the specific heat at con- 
stant volume. Perfect monoatomic gases have 7 = 5/3 
so a convective region (ds/dm = 0) can be modeled with 
a polytrope of index n = 3/2. Using this value in Equa- 
tion (fl]) leads to R/Rq — (M/M ) -1 / 3 and the conclusion 
that fully c onvective stars expand when they lose mass. 

Later on , | Hjellming fc Webbink ( |1987 ) investigated the 
stability of polytropes, condensed polytropes (a poly- 
tropic envelope with a core modeled by a point mass) 
and composite polytropes (an envelope and a core with 
different polytropic indices) for convective (7 = 1 + 1/n) 
and radiative (7 > 1 + 1/n) regions. They evolved 
their models in the adiabatic regime, which means that 
they assumed hydrostatic equilibrium and an adiabatic 
evolution such that the entropy profile remains con- 
stant in Lagrangian coordinates. For the condensed 
polytropes, they showed (their Equation 40) that the 
adiabatic radius-mass exponent, £ a d; asymptotically ap- 
proaches 
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(3) 



where m c is the ratio between the core mass and the to- 
tal mass of the star. Equation ^ describes the behavior 
of Sad m the limit m c —¥ 1 but one also recovers the ap- 
propriate Sad for a complete polytrope for any value of n 
(m c = 0, Equation [T]) . However, the response of stars to 
very high mass loss rates may not be hydrostatic. More- 
over, Equation (p3j) is only valid for condensed polytropes 
which are models that neglect radiation pressure and do 
not reproduce the superadiabatic regime that is encoun- 
t ered in the oute r layers of giants. 



Ge et al. ( 2010 1 also studied the response of mass-losing 
stars in the adiabatic limit but used a detailed equation 
of state instead of a polytropic stratification. The star 
was assumed to stay in hydrostatic equilibrium, and its 
response to mass loss was assumed to be fully adiabatic. 
The entropy and composition profiles were fixed and the 
local value of these profiles during the evolution was ob- 
tained by interpolation from the initial model. For their 
1 M giant star model, the stellar radius increased by 
30%. The study concluded that instability in the mass 
transfer occurs rapidly for donors with a convective en- 
velope, if at all, while donors with a radiative envelope 
may encount er a delayed dynamical instability. |Deloye fc| 
Taam (20101 used a similar approach to study the com- 
mon envelope outcomes of 10 M donors. They found 
that the mass of the remnant can vary by 20% depending 
on when the common envelope phase happens and on the 
initial mass ratio. 

All the models above yield a paradigm in which giants 
expand as a result of mass loss, such that mass loss that 
starts in semi-detached binaries with a giant donor tends 
to be unstable. However, for typical mass loss rates en- 
countered at the onset of a common envelope interaction 
(AT < 1 M /yr, Passy et al. 2012 ), the donor does not 
stay in hydrostatic equilibrium. Moreover, an adiabatic 
evolution assumes that mass loss happens on a timescale 
shorter than the thermal timescale of the mass-losing star 
throughout its interior. 



Recently, Woods & Ivanova (20111 showed that the 
evolution could be locally non-adiabatic, since the 
outer superadiabatic layer of giant stars has a thermal 
timescale so short that it might readjust and reconstruct 
faster than it is stripped away. They present the evolu- 
tionary sequence of a 5 M giant star for various mass 
loss rates, and show that the star grows mildly in radius 
during its evolution (their Figure 3). They also calcu- 
late the critical mass ratio for stable mass transfer for 
different donors. These values are only indicative but 
show that the response of the mass-losing star evolves 
with mass loss and cannot be parametrized using only 
the binary and stellar parameters. 

In this paper we present models of mass-losing stars 
by removing some of the assumptions made in previ- 
ous investigations. We obtain these models using the 
stellar evolution code MESA. The numerical method is 
described in Section [2] We present the simulations in 
Section [3] and verify our method with low-mass zero- 
age mam sequence (ZAMS) models in Section |IJ We 
then study in detail the dynamical response of red gi- 



ant branch (RGB) and asymptotic giant branch (AGB) 
stars, and describe the physical processes involved, in 
Section [5] A summary and conclusions are provided in 
Section (Hj 

2. NUMERICAL METHOD 

MESA is a parallel one-dimensional stellar evolution 
code that uses adaptive mesh refinement and adaptive 
time stepping. In this section, we outline the b asic fea- 
tures of this code. More details can be found in IPaxtonl 



et al. (2011 



In hydrodynamic mode, the full set of differential equa- 
tions of stellar evolution is solved in the Lagrangian de- 
scription: 
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where the mass m is the independent variable and r, 
p, P, T, I, V = dlriT/dlnP, e nuc , e v , c P , V a d = 
(dlnT/dln P) , s, v and G are the radius, the density, 
the pressure, the temperature, the luminosity, the tem- 
perature gradient, the nuclear energy generation rate, the 
neutrino loss rate, the specific heat at constant pressure, 
the adiabatic gradient, the specific entropy, the velocity 
and the gravitational constant, respectively. In addition, 
Xt = (d\nP/d\nT) p and Xp = (d\nP/d\n p) T . To close 
this set of equations, we obtained the equation of state 
from a set of tables computed with the freei^O.^] code, 
developed by Alan Irwin, in the EOS4 configuration. 

In order to improve numerical stability w e use some ar- 
tificial viscosity following the treatment by |Weaver et al. 
( |1978| their Equation 3). Aside from providing better 
stability for the code, artificial viscosity had no effect on 
the evolution based on comparison of sequences carried 
out with or without artificial viscosity. 

Different options have been explored for modeling stel- 
lar mass loss. We first directly set the mass loss rate M to 
a constant value regardless of the evolutionary stage. In 
order to show the robustness of the results, we also study 
how the models respond to a vari able mass loss rate, for 
instance a scaled up wind model ( Reimers||1975 ): 



M = mx4x 10 



-13 
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R 



Me 
M 



[M /yr] 
(9) 

where L is the luminosity and rjn is a dimensionless con- 
stant. We select values of r/R to model very high mass 

4 http:/ / freeeos.sourceforge.net/ 
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loss rates. Typical models for RGB winds use values of 
?7r, around 0.5. 

3. THE SIMULATIONS 

We perform 21 MESA simulations (Table fT). We first 
investigate the behavior of low-mass Z AMS stars (models 
1 to 4) in order to verif y our method th rough a compar- 
ison with the results of 
We then study the ste. 



Ge et al. (2010) for such models, 
iar responseot a 0.89 M RGB 
star (models 5 to 9) in hydrodynamic mode for various 
mass loss rates. Models 10 to 13 are equivalent to mod- 
els 5 to 8 except that they are carried out in hydrostatic 
mode, such that dynamical effects are isolated. Addi- 
tional sequences for a 0.74 M AGB star (models 14 
and 15) and for a 5 M Q RGB star (models 16 to 18) 
allow comparison with the results of |Woods &: Ivanova| 
(2011). We also verify that changing the atmosphere 



boundary conditions from the "simple atmosphere" de- 
fault option to the "Eddington grey" option (model 19) 
does not modify the out come of our simulat ions (these 
options are described in Paxton et al. 2011). In order 



to verify that the initial response of the star is captured 
accurately in our simulations, we finally examine extra 
models (models 20 and 21) which are similar to models 8 
and 17, but with an initial timcstcp smaller by an order of 
magnitude. Time-stepping automatically readjusts and 
we find no difference between the corresponding models. 

4. LOW-MASS ZERO AGE MAIN SEQUENCE STARS 

Low-mass ZAMS stars provide a simple case to com- 
pare our method — solving the ful l set of stellar eq ua- 
tions — with the method used 



Ge et al. (2010) 



assuming an adiabatic evolution (cT. Section [l|) . 
deed, the global thermal timescale of a 0.3 M Z 



star is 



GM 2 
RL 



7.7 x 10 s years, 



In- 

AMS 



(10) 



which is several orders of magnitude longer than the du- 
ration of the different sequences for the mass loss rates 
considered (Table [I]) . Consequently, the star does not 
have enough time to readjust thermally throughout its 
entire interior. Moreover, one can calculate the local 



thermal timesca le of the outer part of the star (Woods 
fc Ivanova|[2rJlT| ): 



tKH,loc(m-) 



M 



u(m )/L(m ) dm 



(11) 



where u is the internal energy p er unit mass and m is 
the mass coordinate. Equation (11) implicitly assumes 



that one can parse the stellar interior at any mass co- 
ordinate into an outer and an inner zone, and that the 
outer zone can thermally readjust independently. This 
might not strictly be the case in a convective region 
where convective overturn is usually shorter than the 
thermal timescale of a given zone, thus stabilizes ther- 
mal perturbations on a timescale shorter than the ther- 
mal timescale of the zone. Figure [l] shows that even the 
outermost 0.01% of the mass needs about 1000 years to 
thermally readjust to the perturbations induced by mass 
loss. Again, this timescale is much longer than the dura- 
tion of the simulations (25 years for the 0.3 M© ZAMS 



star with a mass loss rate of 10~ 2 M /yr, model 1), 
and so the outer layers cannot thermally readjust. The 
evolution of the mass-losing star is thus adiabatic, both 
locally and globally. Among low-mass ZAMS stars, those 
with lower mass have deeper convective envelopes. The 
lowest-mass ZAMS stars (M < 0.3 M Q ) are fully convec- 
tive, and so should be well approximated by a polytropic 
expansion (£ a( j = —1/3, Equation [3]). 
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Fig. 1. — The local thermal timescale (top) and specific entropy 
(bottom) in the outermost 10% of the stellar mass for the Mo = 
0.3 M ZAMS star (dashed black, models 1 and 4) and the M = 
0.89 M RGB star (solid red, models 5 to 8 and 14). 



We carry out simulations where the mass-losing star is 
a ZAMS star with a mass ranging from 0.3 to 0.5 M 
(models 1 to 4). The resolution is approximately 1000 
zones. We reproduce the results from |Ge et al.|2010| quite 
accurately (their Figure 3). For all models, the radius of 
the mass- losing star increases as mass is lost (Figure [2]), 
in particular the 0.3 M model, which is almost fully 
convective and therefore behaves most like a polytrope 
(Equation [I]) . One should emphasize that the small dif- 
ference between the 0.3 M model and the polytropic 
limit arise from deviations of the stellar structure from a 
complete polytrope. We try different mass loss rates of 
10~ 2 and lO^ 1 M /yr for the 0.3 M star (models 1 and 
4, respectively) and find no difference whatsoever. The 
evolution for both mass loss rates is much too rapid to 
allow any local thermal readjustment of the outer parts 
of the star (Figure [I]) . In conclusion, the evolution of 
low-mass ZAM S stars is fully a diabatic. Thus, the ap- 
proximation by Ge et al.] ( |2010 ) is appropriate for these 
low-mass main-sequence stars. 



5. GIANT STARS 

The case of giant stars is somewhat more complicated. 
Giant stars have a convective envelope in which the en- 
tropy profile is flat. For intermediate-mass giant stars, 
the global thermal timescale of the star is still large in 
comparison with its dynamical timescale, and so the stel- 
lar interior cannot thermally readjust during most simu- 
lations. However, intermediate-mass giants also possess 
a cool, low-density outer layer in which convection is very 
inefficient. The local thermal timescale of this superadi- 
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TABLE 1 

The main parameters for the SIMULATIONS: the model number, the main sequence 

MASS OF THE STAR (M M g), THE STELLAR (Mo) AND CORE (Mo) MASSES, RADIUS (-Ro); AND 
LUMINOSITY (Lo) AT THE START OF THE MASS LOSS PHASE, THE MASS LOSS RATE OR REIMERS 
PARAMETER (M), AND WHETHER OR NOT THE MODEL IS CARRIED OUT IN HYDRODYNAMIC 

MODE. 



IV J. (J tic 1 


JlJ MS/ ivl 


Mr. /Mn 




Tin JTt 




M 




i 


0.30 


0.30 




0.28 


1.3 x 10 -2 


-t r\ — ^ TV J 1 

10 M /yr 


No 


2 


0.40 


0.40 




0.35 


2.3 x 10" 2 


10~ 2 M /yr 


No 


3 


0.50 


0.50 




0.45 


4.1 x 10~ 2 


10- 2 M /yr 


No 


■1 


0.30 


0.30 




0.28 


1.3 x 10~ 2 


0.1 M /yr 


No 


5 


1.00 


0.89 


0.41 


102 


1.20 x 10 a 


10" 3 M Q /yr 


Yes 


6 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


lO" 2 M /yr 


Yes 


7 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


0.1 M Q /yr 


Yes 


8 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


1.0 M w /yr 


Yes 


9 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


m = io 6 


Yes 


10 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


10~ 3 M /yr 


No 


11 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


10~ 2 M /yr 


No 


12 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


0.1 M /yr 


No 


1.3 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


1.0 M /yr 


No 


14 


1.00 


0.74 


0.52 


111 


1.37 x 10 3 


0.5 M /yr 


Yes 


15 


1.00 


0.74 


0.52 


111 


1.37 x 10 3 


1.0 M /yr 


Yes 


16 


5.00 


4.99 


0.61 


50 


8.85 x 10 2 


10~ 2 M /yr 


Yes 


17 


5.00 


4.99 


0.61 


50 


8.85 x 10 2 


1.0 M /yr 


Yes 


18 


5.00 


4.99 


0.61 


50 


8.85 x 10 2 


r? R = 2 X 10 7 


Yes 


19 a 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


0.1 M /yr 


Yes 


20 b 


1.00 


0.89 


0.41 


102 


1.20 x 10 3 


1.0 M /yr 


Yes 


21 c 


5.00 


4.99 


0.61 


50 


8.85 x 10 2 


1.0 M /yr 


Yes 



a Similar to model 7 but with different boundary conditions. 

b Similar to model 8 but with an initial timestep ten times smaller. 

c Similar to model 17 but with an initial timestep ten times smaller. 
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Fig. 2. — The evolution of the stellar radius as a function of 
stellar mass in terms of the initial radius (Ro) and mass (Mo) of 
the ZAMS models (1 to 4) with mass loss rates M = 10~ 2 M /yr 
(low) and M = 0.1 M /yr (high). Also shown is the evolution in 
the polytropic limit (Equation [ill. 



abatic layer is very short in comparison with the global 
thermal timescale of the star: the outermost 1% of the 
mass thermally readjusts in approximately one year (Fig- 
ure [lj. Therefore, the layer might have enough time to 
thermally readjust locally, depending on how the local 



thermal timescale of the superadiabatic layer and the 
time needed to strip it away compare. Assuming that 
the layer has a mass TO sri0 n, its local thermal timescale is 
iKH,ioc(w s hcii) and it takes m sri0 ii/-M to remove it. This 
leads to a critical mass loss rate: 



M crit 



iKH,loc("'l s hell) 



(12) 



which gives the threshold for the readjustment of the su- 
peradiabatic layer. If M <C M cr it, the superadiabaticity 
cannot be removed and the outer layer thermally read- 
justs. If M > M C rit, the outer layer does not have time to 
readjust and superadiabaticity is lost progressively. The 
higher the mass loss rate, the sooner the superadiabatic 
layer disappears entirely, after which the star evolves adi- 
abatically. A sim ilar argument has been made byf - 



& Ivanova (2011 1, except that we consider here the time 
required to remove the superadiabatic layer rather than 
the time to strip away the entire star, as they did in their 
Equation 3. There is no unique definition of this outer 
layer, but for our 0.89 M RGB model one can estimate 
from the entropy profile that m s heii is about 10 -3 M 
and tKH,ioc(?7i s heii) ~ 0.04 year, which leads to a critical 
value for the mass loss rate M cri t ~ 2.5 x 10 -2 M Q /yr. 



5 



5.1. The canonical case of a 0.89 Mq red giant branch 

star 

In this section, we first study the canonical case of 
a 0.89 M Q RGB star, and compare ou r results to the 



adiabatic models from Ge et al. (2010). We carry out 
hydrodynamic simulations with constant mass loss rates 
ranging from 10~ 3 to 1 M /yr (models 5 to 8) and one 
model with a varying mass loss rate (model 9). We also 
carry out their hydrostatic counterparts (models 10 to 
13) in order to compare with previously published mod- 
els, and to demonstrate the error a hydrostatic assump- 
tion causes. The resolution for all the one-dimensional 
models discussed here is approximately 2500 zones. 

We plot in Figure [3] the evolution of the stellar mass 
and the mass loss rate for models 5 to 9, while the evo- 
lution of the stellar radius for models 5 to 14 is shown in 
Figure|4j Mass-losing giants barely expand, if at all. The 
difference between the hydrodynamic and the hydrostatic 
models for the lowest mass loss rate (M = 10~ 3 M /yr, 
models 5 and 10) is hardly noticeable. This is due to 
the fact that for such a low mass loss rate, the star is 
barely driven out of hydrostatic equilibrium. A compar- 
ison between the acceleration (a — dv/dt) and the gravi- 
tational acceleration (g = Gm/r 2 ) profiles confirms that 
this model stays in hydrostatic equilibrium (Figure [5]) . 
If one increases the mass loss rate to 10 -2 M /yr, some 
marginal differences arise in the very early phase of evo- 
lution between the hydrodynamic and the hydrostatic 
models. This is the threshold for which hydrodynamic 
effects can no longer be neglected, as the acceleration in 
the outer parts of the giant represents almost 1% of the 
gravitational force at that location. These effects natu- 
rally increase as the mass loss rate increases. 
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Fig. 3. — The evolution of the mass (top) and the mass loss 
rate (bottom) for the 0.89 Mq RGB star with mass loss rates of 
10~ 3 M /yr (dotted magenta, model 5), 10 -2 M Q /yr (solid blue, 
model 6), 0.1 M /yr (dash-dotted red, model 7), 1 MQ/yr (solid 
green, model 8), and variable (dashed black, model 9). 
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Fig. 4. — Top: evolution of the stellar radius as a function of 
stellar mass for the 0.89 M RGB star with mass loss rates of 
10~ 3 M /yr (magenta), 10~ 2 M /yr (blue), 0.1 M Q /yr (red), 
1 Mg/yr (green), and variable (black). Both hydrodynamic (solid) 
and hydrostatic sequences (dashed) are shown. Bottom: a close-up 
of the early evolution of the sequences shown above. 
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Fig. 5. — Ratio between the acceleration and the gravitational 
acceleration after the 0.89 Mq RGB star star has lost about 
10 — 2 Mq. The mass loss rate is (from bottom to top) 10 — 3 M /yr 
(dotted magenta, model 5), 10 — 2 Mg/yr (solid blue, model 6), 
0.1 M /yr (dash-dotted red, model 7), and 1 Mg/yr (solid green, 
model 8). 



Hydrodynamic models for which dynamical effects can- 
not be neglected (models 6 to 9) all contract in the early 
evolutionary phase. The higher the mass loss rate, the 



more the stellar radius decreases (Figure [4]). Later on, 
differences arise as stars with high mass loss rates keep 



G 



contracting while stars with lower M first re-expand 
slightly and then contract again. For the model with 
M = 10 -2 M Q /yr (model 6), the radius of the star 
grows by less than 5 %. We cer t ainly do not see the 30% 
expansion found in Ge et aL] (2010) for a 1 M Q giant 
star (their Figure Of. Model 8, for which the shrink- 
age is most dramatic, considers a typical mass loss rate 
that is encountered during a common envelope evolu- 
tion. On the other hand, the equivalent hydrostatic 
"test" models (models 11 to 13) all show an expansion of 
the radius, with the higher mass loss rates leading to the 
largest increases. This be havior is similar to sequences by 
Woods & Ivanova (120111 their Fi gure 3). This compari- 
son demonstrates that dynamical aspects play a critical 
role in the stellar response of our model. Some energy 
that would be transformed into internal energy or expan- 
sion work in the hydrostatic assumption can now go into 
kinetic energy. 

In order to understand the reasons for these different 
behaviors, we plot in Figure|6]entropy profiles at different 
times for models 6, 7 and 8. The evolution of the entropy 
differs significantly between cases with different mass loss 
rates. For the lowest mass loss rate (M = 10~ 2 M0/yr, 
model 6), the entire interior has enough time to adjust 
thermally. After 40% of the initial stellar mass has been 
lost, the star has still the entropy profile similar to the 
one of a giant star stratification. For the intermediate 
mass loss rate (M = 0.1 M Q /yr, model 7), very lit- 
tle of the outer layer loses its superadiabaticity in the 
early phase. Eventually, the interior of the star has not 
changed except in the outermost parts. Moreover, the 
superadiabatic layer is much less prominent than for a 
regular giant star stratification, and a radiative zone has 
developed just beneath it. In the highest mass loss rate 
case (M = 1 M /yr, model 8), one can see the superadi- 
abaticity being removed very early in the evolution and 
the build-up of a radiative zone. Eventually, the surface 
is not superadiabatic anymore. The entropy profile of the 
interior layers has not changed and the mass for which 
the entropy drops is negligible. Only a small radiative 
zone remains on top of the convective zone (Figure [6j 
middle and bottom panels). 

The evolution of the radius profile as a function of mass 
for the hydrodynamic case shows that while some mate- 
rial close to the surface always moves out, the radius 
continuously decreases (Figure [7]) . A certain kinetic en- 
ergy is associated with this localoutward motion which, 
in the case of hydrostatic models, goes into expansion 
work (potential energy). As a result the latter models 
increase their radius. 

We can try to explain the formation of the radiative 
layer mentioned above by considering a tiny "sub-layer" 
at a constant mass coordinate, located within the su- 
peradiabatic layer of the mass-losing star (Figure [8| . As 
the above layers are removed, the sub-layer can as a re- 
sult more easily radiate some of its energy outwards: the 
temperature decreases while the local density stays al- 
most constant. The thermal timescale is shorter in this 
phase than the dynamical timescale. Later on, the den- 
sity profile readjusts and the density of the sub-layer 
drops significantly while its temperature only decreases 
by a small amount. Now the thermal timescale is longer 
than the local dynamical timescale. These two phases 
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Fig. 6. — Entropy profiles at the onset of mass loss (dashed black), 
during the early evolutionary phase (solid colors, right panels) and 
after the 0.89 Mq RGB star has lost about 40% of its mass (red, 
left panels). The mass loss rate is 10 -2 MQ/yr (model 6, top), 
0.1 M /yr (model 7, middle) and 1 M Q /yr (model 8, bottom). 
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Fig. 7. — Radius profiles at the onset of mass loss (dashed 
black) and during the early evolutionary phase (solid colors) for 
the 0.89 Mq RGB star and a mass loss rate of 1 MQ/yr. Both 
hydrodynamic (top, model 8) and hydrostatic sequences (bottom, 
model 13) are shown. 



are also seen in the various entropy profiles (Figure [9 ) . 
First, the temperature in the sub-layer decreases while 
the density stays almost constant, leading to a decrease 
of the entropy. Then, the density drops while the tem- 
perature only marginally decreases, which leads to an 
increasing entropy. This suggests that the readjustment 
of the star happens in two (nearly) distinct phases: first 
a thermal readjustment during which some of the energy 
of the layers is radiated away, then a dynamical read- 
justment during which the density in the outer layers 
decreases leading to the formation of the radiative layer. 

5.2. Additional models 

We also carry out evolutionary sequences for a 0.74 M Q 
AGB star with two different mass loss rates (models 14 
and 15) and for a 5 M Q RGB star (models 16, 1 7 and 
18) similar to the one used by Woods & Ivanova ( 2011[ ) 
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Fig. 8.— Profiles in the p - T diagram for the 0.89 M RGB 
star (M = 1 Mq /yr, model 8) at the onset of mass transfer (solid 
black) and at different times during the early evolution (colors, the 
total mass is given in the legend). Also plotted is the location of 
a fixed mass coordinate (dashed black with crosses, from left to 
right: 0.885, 0.883 and 0.881 M ). 
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Fig. 9.— Entropy profiles for the 0.89 M RGB star (M = 
1 M /yr, model 8) at the onset of mass transfer (solid black) and 
at different times during the early evolution (same colors as in 
Figure [Si. Also plotted is the location of a fixed mass coordinate 
(dashecTblack, from left to right: 0.883 and 0.885 M ). 



in order to verify that the behavior seen for the 0.89 M Q 
RGB star case is not a special case. The evolution of the 
radius for both stars is shown in Figure [10| and Figure 11 
respectively. Again, all the models initially shrink in ra- 
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dius. Models suffering a higher mass loss rate then shrink 
faster. In the particular case of the 5 M RGB star, 



our re sults differ from the findings by |Woods S Ivanova 
(2011). Indeed, their Figure 3 shows that for all mass 
loss rates except the lowest one (M = 10~ 3 M /yr), 
the mass- losing star slightly expands. The higher the 
mass loss rate, the larger the star expands. This behav- 
ior is quite similar to the one seen in Figure 12] for the 
hydrostatic simulations, so it is possible that the results 
presented in their Figure 3 have been obtained assuming 
hydrostatic equilibrium. Nevertheless, the stellar radius 
after 0.6 M has been lost by the star (see their Fig- 
ure 2) seems quite consistent with ours, although making 
a more detailed comparison is difficult as their Figure 3 
only shows the evolution until the star has lost 0.5% of 
its total mass. Our 5 M Q RGB model has a smaller core 
mass fraction (0.122) than their model (0.171), which 
might also lead to differences. 
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Fig. 10. — Evolution of the stellar radius as a function of stellar 
mass for the 0.74 Mq AGB star with mass loss rates of 0.5 Mq/jt 
(black, model 14) and 1 M0/yr (dash-dotted red, model 15). 



6. SUMMARY AND DISCUSSION 

We used the MESA stellar evolution code to carry out 
one-dimensional hydrostatic and hydrodynamic simula- 
tions of mass-losing stars, either low-mass ZAMS stars 
between 0.3 and 0.5 M Q , or 1 or 5 M (main sequence 
mass) giant stars, with several constant and variable 
mass loss rates up to a few M /yr. 

We hrst tested our numerical method against the l ow- 
mass ZAMS stars case and reproduced the results of Ge 
et al. (2010). Therefore, it is correct to assume that the 
evolution of the star is adiabatic in this specific case. 

We then investigated the case of a 0.89 M RGB star 
for five different mass loss rates. We showed that the 
mass-losing star does not remain in hydrostatic equilib- 
rium for high mass loss rates and that the evolution is not 
adiabatic, as the outer superadiabatic layer has enough 
time to thermally relax. Only for low mass loss rates 
(M < 10 -2 M /yr) have both the outer superadiabatic 
layer and the stellar interior enough time to thermally 
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Fig. 11. — Top: evolution of the stellar radius as a function 
of stellar mass for the 5 Mq RGB star with mass loss rates of 
10" 2 M Q /yr (black, model 16), 1 M /yr (dash-dotted red, model 
17), and variable (dashed blue, model 18). Bottom: a close-up of 
the early evolution. 



readjust. The superadiabatic layer progressively recon- 
structs and survives the entire evolution, making the evo- 
lution not adiabatic both locally and globally. For high 
mass loss rates (M > 0.1 M /yr), the outer part of the 
star progressively loses its superadiabaticity and the in- 
terior does not have enough time to thermally readjust. 
Even though a fraction of the initial superadiabatic layer 
might survive, a larger radiative zone emerges below it 
and the star keeps shrinking during the entire sequence. 
The evolution of the star is locally non-adiabatic and 
hydrodynamic, as some energy that is stored in gravi- 
tational form in the hydrostatic models is actually in a 
kinetic form, leading to the star contracting instead of 
expanding. 

We also carried out additional simulations for a 
0.74 M AGB star with a core mass of 0.52 M Q and 
a 5 M RGB star. These models are consistent with 
our previous findings an d with the 5 M RGB model in 
Woods & Ivanova (2011). We have also verified that the 
outcomes of our simulations depend on neither numerical 
parameters such as the initial timestep adopted, nor on 
boundary conditions. 

According to our stellar evolution models, giants barely 
expand, if at all. This result impacts the condition for the 
onset of the common envelope phase. Using the Edding- 
ton luminosity limit, one can estimate the mass loss rate 
above which a dwarf would be unable to accrete material, 



to be about 10 M /yr. For higher mass loss rates, the 
hydrostatic assumption is violated and there is no ex- 
pansion of the giant's envelope. As a consequence, the 
positive feedback from the mass-losing giant discussed in 
Section [l] may be reduced. Further investigations are re- 
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quired to quantify how this feedback affects the temporal 
evolution of the mass transfer rate. Overall, criteria for 
unstable mass transfer based on adiabatic mass loss mod- 
els should be re-investigated. Moreover, if giant stars do 
not expand as a result of mass loss, this process does also 
not help the envelope ejection during a common enve lope 
interaction, as speculated by De Marco et al. (2011). 
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